


Institutional Archive of the Naval Postgraduate School 





Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1970 


Energy loss of high energy electrons tn tin, 
lead, and gadolinium 


Mosbrooker, Michael Lee; Sandquist, David Lee 


Monterey, California ; Naval Postgraduate School 


http://ndl.handle.net/10945/14987 


Downloaded from NPS Archive: Calhoun 


Calhoun is the Naval Postgraduate School's public access digital repository for 


f (8 D U DLEY research materials and institutional publications created by the NPS community. 
«ist : Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 


NY KNOX appointed — and published -- scholarly author. 

| | LIBRARY Dudley Knox Library / Naval Postgraduate School 

411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http://www.nps.edu/library 





- 


ENERGY LOSS OF HIGH ENERGY ELECTRONS 
IN TIN, LEAD, AND GADOLINIUM 


by 


Michael Lee Mosbrooker 








United States 
Naval Postgraduate School 


ENERGY LOSS OF HIGH ENERGY ELECTRONS 
IN TIN, LEAD, AND GADOLINIUM 


by 


Michael Lee Mosbrooker 


and 


David Lee Sandquist 





June 1970 


This document has been approved for public re- 
Lease and sale; its distribution is unlinited. 


1136107 





Energy Loss of High Energy Electrons 


in Tin, Lead, and Gadolinium 


by 


Michael Lee Mosbrooker 
Major, United States Army 
B.S., California State Polytechnic College, 1960 


and 


David Lee Sandquist 
Major, United States Army 
B.S., Iowa State University, 196] 


Submitted in partial fulfillment of the 
requirements for the degree of 


MASTER OF SCIENCE IN PHYSICS 


from the 


NAVAL POSTGRADUATE SCHOOL 
June 1970 


ABSTRACT 


The LINAC at the Naval Postgraduate School, Monterey, was used 
to accelerate electrons to energies ranging from 52 to 92 MeV in 
order to study the energy distributions of high energy electrons 
before and after passing through layers of tin, gadolinium, and 
lead. The thickness of these materials ranged from 0.8 to 5.9 g/cm*. 

The most probable energy losses agreed with the theory of Blunck 
and Westphal for all materials used, while distribution half-widths 
agreed only for absorbers of thickness less than 3.0 g/em?. The 
thickness at which theory and experiment began to exhibit a notice- 
able discrepancy was found to be dependent on the atomic number of 
the material. 

Where comparison was possible, results of this experiment generally 


agreed with the findings of similar works concluded previously. 
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IT. INTRODUCTION 


The theoretical energy distribution for an initially monoenergetic 
electron beam which has passed through an absorbing material was de- 
termined by Blunck and Westphal [1]. Their work modified previous 
theoretical calculations of Landau [2], Eyges [3], Bethe and Heitler 
[4], and Blunck and Leisegang [5] by adding radiation losses to the 
previous predictions of ionization losses for electrons when passing 
through thin absorbing layers. The distribution assumes energy losses 
small compared to the incident beam energy. The theoretical treatment 
is presented in Section II. 

Several measurements of energy losses by high energy electrons 
passing through various materials have been performed in the 10 MeV 
to 150 MeV range. Of note are the works of Breuer on aluminium [6]; 
Bumiller, Buskirk, Dyer, and Miller on aluminum [7,8]; Goodwin on 
copper [9]; and DeLeuil and Raynis on aluminum, copper, and lead [10]. 
Except for Breuer's work on the Darmstadt linear accelerator, the 
experiments were performed on the LINAC at the Naval Postgraduate School. 
These experimental results are in general agreement with theory for thin 
absorbers (s 2 g/cm“), but discrepancies with theoretical predictions 
are reported for the thicker absorbers. Such discrepancies are expected 
for the thicker absorbers since the energy loss is no longer smal] 
compared to the initial beam energy. 

In this thesis the energy loss measurements are extended to include 


tin (Z = 50) and gadolinium (Z = 64). These metals were chosen to 
provide experimental results for materials with atomic numbers between 


copper (Z = 29) and lead (Z = 82). It was thought that analysis of 
these materials would provide more conclusive information on the 
effects of atomic number on the agreement between experimental and 
theoretical half-widths as reported by DeLeuil and Raynis [10]. 
Thicknesses used ranged from 1.485 to 5.940 g/cm? for tin and from 
0.814 to 4.83] g/cm? for gadolinium; nominal beam energies were 52, 
75, and 92 MeV. Thick lead absorbers (2.825 and 4.236 g/cm) were 
investigated as an extension of the work of DeLeuil and Raynis, and, 
in addition, their thickest aluminum absorber (5.574 g/cm?) was 
investigated as a continuity check between the two experiments. 

The energy loss distributions are characterized by the most probable 
energy loss and the half-width. The half-width is the full width of 
the distribution curve at half maximum. These quantities, obtained 
from experiment and theory, are the basis for comparison between ex- 


perimental and theoretical results. 


II. THEORETICAL CONSIDERATIONS 


The Blunck and Westphal theory of the distribution for the energy 
loss of a beam of monoenergetic electrons in passing through a layer 
of absorbing material assumes that the energy loss Q is small compared 
to the initial beam energy, E.. Let W(Q)dQ be the probability of 
energy loss between Q and Q + dQ, and X be that portion of the loss 
Q due to radiation. Hence, the ionization loss is Q - X. Considering 


these two loss processes, the probability of energy loss is 
Q 
W(Q)dQ = F Wy (Q-X)We(X)axdQ (1) 
0 


where Ws and We are the energy loss distributions for ionization and 
radiation respectively. 

The Landau equation[ 2], as modified by Blunck and Leisegang, is 
used for the energy loss distribution due to ionization. The distri- 
bution is expressed as a function of Landau's dimensionless parameter 


h as follows: 


Can (x - ay)? 
W, (Q)dQ = o(A) dA = ¢& exp [- Se dPy AZ) 
nlf Tey Dia a 
n [ n 
ae E. 
evel, eee pee = 1.1 (3) 


aR aR 
The terms used in equations (2) and (3) are defined as follows: 


Cis ee and 4, are constants given by Blunck and Leisegang [5] and 


are used to fit Landau's distribution to a sum of gaussian functions. 


R is the absorber thickness in cm. 


The quantity "a" is a function of the atomic ntimber Z, the atomic 
weight A, and the density o of the absorber; and g(=v/c) of the 


electrons: 


a = 2 bee MeV/cm. (4) 
A 8 


The quantity b@ is a correction to Landau's theory given by Blunck 


and Leisegang [5]: 


2 3.0 Lan 2 E. 
Dl een 
(ene ) 





ee, 2 5 (oy) 


r 
aR ” 


where the summation is over the ionization potentials of the atomic 
electrons, and ie is the number of electrons with ionization potential 
[i 

Q is the average energy loss due to ionization (no radiation) for 


electrons of incident energy E., and is given by Sternheimer 2. 


and 13] as follows: 


ea 
Q=--, |B+0.43+ In E, - 6 
B 


2 


INO 


m 
-C- a (X, ” 10g, gp/me ) MeV (6) 


where t is the thickness in g/cm, and the constants A.» Basie Xs aes 
and m. are parameters of the absorber material. These parameters for 
tin, lead, and various other material are listed in reference [13]. The 
parameters for gadolinium are not listed but were determined by the 
following method. 

A. for gadolinium was obtained by extrapolating from a plot of A. 
vs. Z/A for the absorber materials listed in reference 13. The result 
was - = 0.0624. B for gadolinium was found in a similar manner from a 
B vs. 1¢ (ionization potential) plot where I was determined from a Z vs. 


I graph. The result was B = 13.3. C was determined from a semi-log plot 


14 


of 1*A/pZ vs. -C, and the result was C = -6.80. The values for aes 
Xy: and m. Were more difficult to determine, but fortunately they 
contribute little in determining Q. These constants were numerically 
determined. The values for the term a. (X)-109 9(p/me)) ® were 
calculated for tin (Z = 50) and tungsten (Z = 74) in the p/mc range 
(p is momentum, and m is the rest mass of the electron) of the energies 
used in the experiment. Using Z as a basis of interpolation, the cor- 
responding values for gadolinium (Z = 64) were estimated for two 
representative p/mc values (100 and 150). The value for X. was 
arbitrarily set to be 3.0, and the a. and Me values were determined 
from a simultaneous algebraic solution based on the interpolated values. 
The results were a. = 0.418 and mo = 2.10. The error in the determina- 
tion of Q using these values was estimated to be 1% or less. 

For Wes the energy loss distribution due to radiation, Blunck and 


Westphal [1] give: 





> aR d 
Ws (Q)dQ = BaR(Q/E;)** (7) 
where 
a = 1.4 x on 2 43 In 73 + v3 cm”! (8) 


and B is a normalizing factor = “GREAT : 
The distribution of total energy loss according to Blunck and 


Westphal is obtained by putting equations (7) and (2) into equation 
(1) and performing the required integration. The result is the energy 
loss distribution for a single electron of incident energy E.. For 


comparison of theoretical and experimental values, this distribution 


function was used, with corrections to account for the finite energy 
width of the incident electron beam. This treatme:t ts descriind in 


Section IV and Appendix C. 
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TTI. EXPERIMENTAL PROCEDURE 


The LINAC of the Naval Postgraduate School was used to obtain 
electron beam energies from 52 MeV to 92 MeV. This beam was elastically 
scattered at 90° from a thin (3.5 mil) aluminum scattering foil, passed 
through the absorber, and finally analyzed by a 120° magnetic spectro- 
meter described by Kenaston, Luke, and Sones [14]. 

This general experimental arrangement was similar to that used by 
Miller [7,8] and DeLeuil and Raynis [10] with the exception of the 
removal of a 3.5 mil aluminum window and the installation of a ten 
channel coincidence counting system. These changes are discussed below. 

The aluminum window removed was located just before the scattering 
chamber. Since the presence of the window caused a broadening of the 
incident electron beam distribution, its removal decreased the energy 
width of the incident beam for this experiment, as compared to previous 
works done at the NPS. 

The ten channel counting system consists of ten front counters and 
a single backing counter operated in coincidence with the front counters. 
The entire ten channel system has an energy spread of about 3%. 

The absorbers were positioned approximately 3 cm. from the scattering 
foil as recommended by DeLeuil and Raynis [10]. Absorber thicknesses for 
tin were 1.485, 2.970, 4.455, and 5.940 g/cm? : thicknesses for gadolinium 
were 0.814, 1.610, 3.221, 4.026, and 4.831 g/cm? : absorber thicknesses 
for lead were 2.825 and 4.236 g/cm. 

A thick aluminum absorber (5.574 g/cm) was used to correlate the 


results of this experiment to the results of DeLeuil and Raynis [10], 


7 


which were obtained before the ten channel countin, systen wa. .: called. 
The results agreed within experimental accuracy. 

Since the electron beam scattered into the absorbing tiai.: was 
not monoenergetic, the energy distribution was determined both before 
and after passing through the absorber. Various thicknesses of absorber 
material, including no absorber, were positioned on a remotely controlled 
ladder device. For measurement of the distribution before passing through 
the absorbers (called a zero peak), the ladder was positioned such that 
the beam passed only through the scattering foil. After measuring the 
zero peak spectrum, the various absorbers were moved into the beam to 
determine the energy distribution of the beam after passing through each 
absorber. 

In the case of tin and lead all absorber thicknesses could be measured 
without turning off the accelerator and thus possibly altering the char- 
acter of the beam. This was not possible for gadolinium. Because of its 
cost, only a small quantity was purchased. Hence, only two gadolinium 
absorbers could be run without turning off the accelerator and rearrang- 
ing the gadolinium on the ladder. Since the beam character could change, 
a zero peak measurement was made whenever absorbers were replaced. 

The data represented the number of electrons detected by the coin- 
cidence counting system at the exit of the magnetic spectrometer. A down- 
stream Secondary Emission Monitor was used as a standard for normalization 
purposes, in that each data point corresponded to a given integration 
current, i.e., a certain number of electrons passing through the absorb- 


ing material. 
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IV. TREATMENT OF DATA 


The data reduction for the ten channel counting system is standard 
for multi-channel systems and is on file in the NPS LINAC computer 
library. In this reduction process, three characteristics of the 
counting system are required. The characteristics are: 

(1) the energy spread of the counting system, 

(2) the energy seen by each front counter, 

(3) the relative efficiencies of the front counters. 

With this information the actual energy corresponding to the front counter 
data is calculated. Counting rate and background corrections are also 
performed during this computerized process. 

In the lower energy regions of the spectrum, where all counter read- 
ings should be nearly equal, a counting discrepancy was noted in that the 
lower energy front counters read consistently higher than those on the 
higher energy portion of the counting spectrum. This effect occurred at 
energies lower than required to obtain distribution half-widths and thus 
did not affect the data reported. However, it made confirmation of the 
shape of the distributions at very low energy (< 20 MeV) impossible. 

The electron beam energy is not monoenergetic as required by Blunck 
and Westphal theory. To compare experimental and theoretical results, 
the zero peak data must be used to modify the eheomasieal predictions. 
With no data treatment, the theory will predict half-widths which are 
too small and, in the case of asymmetric zero peaks, erroneous most 
probable energy losses. Removal of the aluminum window noted in Section 


III narrowed the zero peak to the extent that it could be treated as a 


ie 


symmetrical gaussian distribution. The IBM 360 compiiter of the NPS 


was used to unfold the zero peak energy distribution into the theory. 
A histogram method described in Appendix C was used for this untolding. 
The computer program used to accomplish the unfolding is Appendix D. 

In the comparison of experiment with theory, the measured incident 
energy distribution was unfolded as previously described, properly 
normalized to the experimental data, and plotted. From these plots 
the theoretical half-widths and most probable energy losses were 
determined. The experimental data were plotted along with the correspond- 
ing theoretical curve, and measurable parameters were compared. The curves 


are shown in figures |] through 38 of Appendix B. 
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V. RESULTS AND OBSERVATIONS 





The values predicted by theory and the experimental results for the 
most probable energy loss and half-width for tin, gadolinium, and lead 
are shown in Tables I, II, and III respectively. In these tables, BW 
refers to Blunck and Westphal theoretical predictions, a represents 
the most probable energy loss, and HW is the half-width. Data are given 
for predictions with and without beam folding. 

On comparing the theoretical predictions with beam folding and the 
experimental results, it is seen that good agreement exists for target 
thicknesses less than about 3.0 g/cm for both half-width and most 
probable energy loss. At thicknesses greater than this, the agreement 
in most probable energy loss is still reasonably good, but there is poor 
agreement in half-width. The greatest variation in most probable energy 
loss is 8% while the average variation is 3%. These percentages, while 
averaged over all material, are typical and not material dependent. The 
average half-width variation is 9% where half-widths were obtained. As 
can be seen from Appendix B, figures 2, 4, and 6, the theory does not 
predict a half-width for a certain target thickness for each element. 
This cutoff thickness is smaller as the atomic number increases. For 
example, the theory would predict a half-width for 5.6 g/cm® aluminum, 
but not for 4.2 g/em lead. 

It is concluded from these results that the Blunck and Westphal theory 
predicts correctly the most probable energy loss for all targets used and 


predicts satisfactorily the half-widths for the thinner targets. Thin 


21 


targets can be defined empirically as those satisfying the relation 
T(z) 1/8 < 13.5, where T is the absorber thickness (g/cm?) and Z is 


the absorber atomic number. 


22 





+ |+ |+ |x 


+ | |+ J+ | 


+ |-+ |-+ 1+] 


| ZUusWLUadxy 


+ |+ J+ J+] 


+ |+ |+ J+] 





pap|o4 pap ,o4 
Md Ma [e}uewLuadxy Mg 


(NOW) MH Ca 






UL} JO SOLZSLYazZIeUeY) UOLINGLUISLq $so7 ABuauUy 


SATaVL - W XIONAddY 


pauLezqo 7OU YAPLM 4LeH 


tet Saly 


Md 


‘I FTV 





06'S 
SS? 
Oke A 
S87° L 


wo /6) 4 


0° v6 


Oa 


CS ds 


(A@W) 


~™ 


3 | 


as 


‘a_Bue Yobuey UL UOUUS 48 02 ANP pasn ssouyDLu SALI9I9IIFI xx 


pauLezqo OU UIPLA SLPH x 





RE ape. 


0°02+ 00°S x * 
O°'2 + O'Cl x * 
fp Vel S O53 Lv’s 
ee ee | 96° 78° 
80° + 79° 09" 9¢° 08° 16 
00°SZ | 
Sl os 
08°¢S 
PSP LO Pep |O4 \ 
| e}UuswWLuadxy Mg | eJuswWLUSdXy Mg Mg (_w9/5)3 (A®W) 3 
(APN) (AeW) “d 


WNLULLOPeH JO SILFSLUSTIeULYD) UOLINGLAISLG sSso7 KBusu3 “LI 31gavl 


24 


‘3|Bue z9bue2 UL AOUUD 8 0} BNP PASN SSOUADLYZ SALZIIIIF xy 


paulLezqo 10U YIPIM SLeH x 





OL° + 8°8 
O'S + OLE 
OL" + 8b°9 
O'S + 9°22 
OL" + 09°9 
Pep oJ 
| PJUdWLUadxy | eJuaWLuadx3 Mg 





a a em ea 


(A2W) MH (AQW) 0 


- ve - ~ oe ieee. ee a ee ere ee on 
ee - 
om 
SS RE SE I RA IR LT TE GL SP PS ae, I-A GR ad es ee eS EA pe ee ee - ee — 


pea] JO SOLASLUaZIeURYY UOLINGLuystq sso7 ABuaUy “TIT AMV 


Z5 





Ob"SY U 


Figures 


APPENDIX B. 


is Se ea ear eer errr ery yyy Pian 


Co 
LD 


| 


0 


OQ 


J 


'D 


fo bake} 


z WOW 
OUD 


[| aunbL4 


SSS Oat. elec 
Brae OUP *¢ Se U8 “0 


SAAW co=l4 
iWaxa ° ——  G@2URis Vote ace es 


San TEI aaSsebeee aia 
"SN SSO] ADGANA 4 idbadOdd 1LSOW 





UO °0 


Cae 


a) 


QO "6 QO 
dD 


wie | 


QO°ST 


9 
u 


AAW 


26 





[Pe veatee 


TTT Se eee iy rit tq 


i—*— 


cHWOZW9 
O9°S OB *p U0 *P 


/ 
/ 


Vibe Ns tis fi eae adit aes ba holon alual |} 


Z:a4nBbL4 


ge SNe i caer sales lies ay elt 


Bree 2 OP? Oral UE ‘OU Uv *O 


pot cone rs PTT TTT yet Se 


Leia” == Ga Me ee ee 
‘lhaad° ——* a0 4 oe | ee Ses ie 


hdad =~ WeORt Se erie sree 
SSSINCANMEL. loishsichkisista, (si Ik 
SCs ee 0) ieee 





o0°e2 Oo'T2 OO-bT HO 
ME 


wy Se 


ASW 


2) 





OD"S 


eee era ae eer er rrr rye 


~” 


€ aunbly 


2WI“WO *«SSANHQIHL Ladddl 


DU "PD Tee ODS Oo 


4! 


“i 
I 


et 
SS 


Wdx 4° — 040 104 Me 
Wdxd*—— O40 104 Md 
Nes ———  CeGaeimnic 


UE “0 


lah vee ies 
el eco 
Ne ae 2 Sees 


ANWSTHL dsdad0Sdb WAINT I2da9 
"SN = SSO1 ASHANA 43 ldcbdddd LSOW 





Ue 


UU "0 


CGie 


Us. 
dU 
28 


Bua | O0°& 


aS | 


NAW 





OF *Y 


ee 


U8"? 


WIA 
OU'P 


py ounbL4 


OSS itel ea ee mee elneitalat 


Seema eesti 





OC *E Ue O9'T OB *D Ov'D 
TTT TT 

"I 
j* . 
| 5 
/ 9 
a) 
INJ 
- 
Cy 
r=) 

'tWorxs' —-'OS0 104 4H snsah co=14 
‘tHok4° —:- O40 104 WY fASW Sc=t3 ne 
‘tadeks’ —'OS0 104 WY HANSA cSt) cs 
SSSANMUTHL = YS&HOSEU WONT 180bU = 


wo Se 


MH 


s 
a 


Ash 


29 





OP "Ss 


2 WO“WS 
OPO 2 


4 
et ‘ 


G ounbl4 


MSS aN ei Eee icaie 


OC "E UD *<¢ Gye) DS *O 


| { 


—_ 

— 
~ 
~ 


1 ook 6 eee 1 { 
Peper NCU OE ech ris tO c Gta Ucar Octet es ' es Sa nee ' 


ay . 
& 









¢ - iWGxs” ——=-C200DS Wee eeae 
©: lWaxi* == GS0004 4S “anne Za ke 
*: | Woxs” = G2080e oS) Seema 
SSANWOTHL = dsedbS&s Cbs | 
"SN =SSO1 AQdSNA J leddOdd 1SUW 





OU°O 


Dee? 


QU°S SUrse 
‘dv 


OU "6 


Be Si 


JIN 


= 
_ 


) 


30 





as 


Sa 


UR *P 


cWHJ-Wo 
UU '*P 


SS Se eects melee eae 


9 aunbLy 


Oe 2 


Pane 


ee @ «© «+ # @ 


C2 


QO °2 = 


DUO*DT 
Mei 


See 12 


as 
a 


Ash 


3] 


[——*¥—$— 
ay 


t!iWd¥3i'—:090104 Me ‘Nah c6=i4 
J | ®t iWak3'—:090104 MG ‘Nd Sd=14 


j head: == UAU iS Sere 
GS SNA tee ee ote te ts ta 
oS emcee) rlet atte ig 





i 


lethal festa ia Pies ta Lata EI ie Algae Ada eli lief alii ee | 


OO e2 


St 





aie ists 


Z aunBbl4 





Sees ae 
ci towd 


I 





(ah la éo-la WIZWS GSopt dts 


a Soto Tap 


Bet pti yy ee 


32 





cs 
Ee 
I 
“I 
an 


Hoe? Ce CU at US CP COOP 














Bivintal@in 


ay 
Goll 


aNOGd dd 


2 


So 





2 a 

TSiNSHEdGY 4 ——‘AbGSH 
HIGHUSSH N 

NOT LAGEXLS LG 








6 aunbL4 


ah 
UU PS Ju US JU "SP UU" CP UU "YE UU "Pe. 


i 





ee eo alae 


Fatale deal sbi nile banal 1 | Nee es Pete det! Cae 


UD “Ur. 


AW SS a 

A NOTES) AUST a al smear ei Pern 
oa eed ett 

NOD LOYDELS to 


eo} 
3 
wa 
Cx 


aAdd4nd 


ligwWON 


to 
i: 


SIND Odé 


fotos) elie eee 


34 





OL eunbty 


ah 


UD "SE UO "8d CG 


oa 
i Se 
i 






eae Sale 





d4g00SSb Nid 
NOCLAGTHLStQG xvadydna 








$ @#@ #4 @ «©  # f 


WIS OPB*GiL 
4 HLHO TWiNGUEHdW¥3 ———:AHO3HL 





Sahel 0421 bhdON 


Niel aes ake a es beie tats oleae ot ot 1 fi I Vathee ho Oe Se St een 


35 





Es 


lala til 


bo 


ce 


We 


od 


~ 


C 








LL aunbiy 


rau 


ee Cee tae OiSer0ra Pree 








Nah Gb'peld dbo Geet bes 

Ata) ASI Se el Gli eee lel at 
H4dbubSao NIL 

ei) Sl S| Soe sist 


ee a a 2 — a eR 


| bol ee Pee lt tee lcs 





tad OM 


EeNiea Bazi 


36 








ee ee 


* 


te ke a 


Laciaiteal 


| 


etilicshe tl 





I 


1 i NIN S6*Pd 103 





zt aunbly 


“El 
Oboe OG'6S 


et LOG IWiNaWiddue J 


6 ea ewe ae a 


A Gleneeena 


Se 


HdGuoses NTL 


NOT LIAYTHLS LO 


ti yi ee ell 


ay aN 3 


ell haphlislear Sieh sl ee caine lea 


O42? Wohl 


UD SY 


LN 


eS 
> alt 


37 


an 





’ 
so Ve ee Te ee Oe le A OER US CO 


peed 


BS Sa a a aa a cat ee a elevate dati ee ll 


Sig as 05°03 Uae aes 


€L eunbe4 


eS Seo ve Oe ee ee eee ed Ra fe 


Nan Gb*pail4 IZM SSP PE 
A VCS St x S cceererereee ar el eet 


d4gdbSdo Nit 





NOTLAYTHLSTQ adydN4 


| | | tele eal eel EUS Teta) 


SiNNOD 9421 TbhwdeN 


38 





+ 2S 2s 


(ealtaatinal 


{ 


{ 


Cee: 


[+4 


ppt 


00 "SS 


pl eunbLy 


4h 


OU "PP Oe e 


NdW Sb*pa: td 


A YLGQ ItiNaWlusdas 


Homes oes U 


ef i ee we we GS 


Telia siete 


Sa eel dl 


YANGUSHE NCL 


NOL LNE TELS CG 


AJHINS 





421 TbWdaNn 


SINNOD O 


39 








G 


6 


GL eunbey 


m4 
OS‘T6  o0°TE 
1 +? 
I 





rs 


CE “bot [4 
aOibG ldaNaWlusae 3 





Gee 


ee 
ey 
oA 


eae ee 


cb 


YaubbSuo Nil 


NOC LINE LeLS tu 





asd 


UU * OG 











es 


ao 


ON 


Toh 


SaNAOl dg 


40 





bo ee ee te a 


ee eee ai eer rey 


| 


Bhatia 


UU "U6 


LIA 


caane. 


g{ aunbly 
al 
DS 63 Uae C 


' i 
} ‘ ‘ ‘ 4 : ‘ 


Gene Ce as 
Argiod IhiNaWiddaea 


d4GUD054b 
NDT LA YlSstu 





Wi-Wo Uch ct 
eae ape 


NUL 
ADEs 


gles ealied holt ee! 


C2 


TiohidiN 


ai 


a 


GSaNPe? da 


4) 





Oe 


Sr re en ee Pe | 


De eee) 


ae 
ak 


Dias Jia oti a Miat a Uiced ca ea 


U0° 88 


+r 


t--- 


O0 "PS 


LI 


00°08 


i 
WEL 


aunb 4 


sh 


UO Se ee 


c& "Pb tts 
GU IbLNAW odds 
HIGHS NIL 

NOT INE CELS tO 


| B22 ee 


DO°8Y 


‘ é ‘ é i 


c(HI-WI GSP PEL 
ee SE eee 


yale ee 





aS 


KIN 


| 
See! CAA 


42 





Cee 


Z SST ee TTT TTT Ty 


{y 


Li 


rcs 


DU "PS 


BL aunbL4 


NAW 


OU" 9S UU" 2D UD Yd 


Ada co perls 
X°OLIE ISLNaWldddax4 
didaboads NIL 

NOTING THIS IG 


Ze avis 6 ioe eieee 


eee el eta 


AJHANA 





So 
C 


43 





eS 
Cy 


eee TPP UCTTereey Terrier 00 


- 


4 





—] 


6l aanbly 


AAW 
DUT TS Os, as OU “OS 


eee oo le 


Kolb IbiNsaWldsdx4 


H44d0So6 WIND IGb5 


TOUTE STE apenas) ee 


1° GP 


2WI“WS DIS" OeL 
"Ad SHL 





CU 6D 


| IbWolN 


4 


Sm) etlaZ 


44 





E 


halal 


6 6 «@ #¢ @ 8 4 


ae aoa 


islets talagedl: | 


lait! 


02 aunbty 


na 
aaa Oeta2 UU "bP US "ep 


ia 


Lae ie 


Eo eae Saale ae tee eaten 
AU GS Siete els ocr ee ele ate 


d34uDS6o WAIN I@ubd 
Pi IA Veta teres aS esnips ae 





5 » a 


SIANNOD O42é1 lbhdln 


45 





ee 


(Aled a [algae 


{ 


eos 


eA x ar 


! 


hat 


0G 8p 


U0 "dp 


ef ce 8 SU Ce 


[2 aunbly 


rd 


OB ea? U0" PP ere 


Man 8c"dS: [4 


A‘WG LAU IWiNAWC ddd 


cH WS 


AdL 4! 


HIGYDSEH WENT Ta0HS | 
NOLLAGIULSIG AQYSNI 
| | } 


va I Mie 


¢ 






aah 


SINNOD OAél TbkYON 


46 





<< 
fe 
iC 
ees, 


™“ 


| 


| 


go aunbl4 
4 


UP UU "OP ee 


Seti etey ene 


UL "Pe 


4 4GHOS 4b 


IW GCUCRs 4 


ea el 1) = | etIE ier oe HG ata ered a cee aE 


WATN( QUE 


| IbROEN 


7 
i 


47 


Sai) iidé 


NOL LANETYLSEG 


ATHIN 4 





fe td Pe ee at eae 





bay able aia 


¢ © t+ «@ mw 4 


Vee 


or | 
rj 


DU*OP 





€z aunbly 
sh 


Cl le DU" Pb? DUST UU"? U0 "0 


& € A €  elUmeDhUCK!UCUGUlUlU er ies ee 


flea valaees! letersi 


dW) Tee "Pel 
Aaa US Sale eee eee ap 
HAGHOSHE = WATNT 1OU89 
a ee eC SRE cL Seat Caras ast et 


SINNOD O4é1 IbhdoNn 


48 





ye aunbL4 


Ns 
Ce OS ae eG Uae Vee a ee Ue Ce OS Ne OU'Te¢ 


» CECT TE el a ee see ee ae al 


ae 





SINNAO2D 0421 TWhWyON 





(an Sere eWJ7ZWO PIB Ot 
oy SC) ASIST RSIS ale 2 eee rane een 
HAGHOSEHE = WATNT 0082 

NOTLAGTELSIG xvQYsNS 


ee ee ee ee eee 


Pewee eg dette 


49 





UU “pe. 


Die G 


tI 


— 


Pee aly 


JU 


ue 


oa 


| 


US 


3 


G 


a 


G2 eunbt4 


Aah 


ae OG ie UO * i 


aaa 
faa 
fo] 
t= 


Nah Yo*pa'l4 
AtOLYG TWiNGWiddaa 3 


H4eYOSHb WAL 
NOLTLAYTYLS LO 


WIZ OTST Hh 


ey pe aia 
NO 1308) 
ASHAN 4 





\ 
J 


iaWdh 


i 


SINMOD dA¢ 


90 





gz aunbl4 


dn 


Vee Uae OS le OIG) es areas 00°39 Cee DES ga 


eG ee we 


(a eae ta) Ales inl  Care ge 
SU Ve Se elle (eles ee Oe 
U4GHDSH = WAENE IbUbo 
NOELAGTHLSIG aQd4Nd 





sll eisalit Ie ae SSE ct a Fe eee { | l Le a aie ates a inal) a ena 


SAiNMGD JAZ] Wohin 


o1 





J0°88 


| 


Cac 


bk ee 


beg 1 
1 ed 
! st 


Oa 


| 
HFT 
th 


ty Ki 


aig I 
ie 


= 


eat | ee 


Opes 


Z2 eanbly 
4h 


U0" PP OD ele OO SeE Oly SL 0D 'D 


i 8 Ere be Wy ea aes Cn. cea ee cares 


4 


~ 


NdW Soha: fs aWI7W5 SCOR EL 
eC) a TS el ee eres tose 


HAGuUSH WALNE (Odd 
ROSES eee 





SINNOD 0421 IoweONn 


02 





UU*GE ee at 


leslesl alee let 


Hes 


aed 


zah bets I 
HE 
? 


inne 


[het 


[este cabal es peal 


i 


Wa 


— hes 


atte 


ces 


8z aunbly 


Nah 
Se U0" PP Ce OU 


et 
}-~-- or 
tbs 1 


M4 Yo"re ls 
A POLE In iNdWl oda 4 


EG OO"TI 


eWIZW Tee ped 


—— ree 


UAgHOSHR WAIN G00) 


NOT LAL ELS TO 


alee tes 





' 4 


SAaNN@2 C42éi TbheON 


23 





6z aunb: 4 


ah 
UU "Tb Serie OU “Ob Oe 00 "be Cee ae UD ' 8d 05 "eu 


t+ 
= 
4 
Se 
4 
! 


+ 
tH 


: Pe im areola 27 essen 
i PAA VST Staal Sere TU 


Y4yuebSd6 WANT 1O0bS 
NOTINGTHLSitO aAgdaNns 


Heels UNDE beac Ticad Debaebetc 11] 1 jl | Ret 1 MST edie i Nai ieee 





NAB C421 IbheGN 


i 


S 


54 





Cio 


O€ sunbly 


Mahi 
Peete 00 *Tb 05°06 00°06 Opes DO ‘b8 HepeeS 


# 
x 


a 


Naw be lets AWIZWO PIB CE! 


LO Ve a sl Sil te rc eh eee 


dd6yOSEa = WATNE Wd 
Pieletet NSD SG IN) ie leet 





~ 


42] IbWYWON 


SINMODD dG 


99 





L€ aunbLy 


mh 
Pee U0 "68 Seco 00°98 UG" PY Ota Pee OU "OR oe 


Qe SE OY ear ee ee 


Ndi pe*tor ts EA 2 lt) eae oc 

a M6 YLg0 WLNSW[ esdk4d ———‘ald4aHi 
< H48u0SEG WENT duo 
SCE NS eS) Sera eletae 





Wiad bs Deca peal mk a tint ce cae ee Nee in| | | |Jobs ok oleae! 


SINNOD d3é1bkHON 
96 





ze aunbly 


M4 
UI Gi OU °8b OURS Oe HU SIS O0'CP 00 °8¢ OU'PT 
: i 
i 

AT 

1 
ay, 
hi Nah Pe lo:[4 ea 7 ali delle ei 
F A'OLGU WiNshW [dda i —m vnddahi 


YAGYDSHE = WATNE Ob 


~ NOTLAYLYHISTY 


AddsN4 





UU'U 


SAiNNOD d3é1 lendONn 


2/ 





COCs 00°86 OO ° PS 


eet eo ald! TTT | ala et 





C0) KG 


\ 


€€ aunbl4 


Nasu 
00°SS Oe ae 


00 °8¢ 


so. @# # 8 @ 


pei 
a 
HEY 


Pe ZGueee eee WIZ C8e "het 


cL Ve UNE Clete cnet eee aed 
W4gGyOSEd WINE :Od02 


NOCLAGLYLStG aAdYsAN4 





ere io 


SiNNO9D 0321 TibhdeNn 


08 





Gs Oe 


a7) tee a Sees 38 


ro 1 ttt 


US 


ae 
rs 


DS oP 


O00 °8P 


ye aunbl4 


nau 
OS'9 00'Sh OS'Eb 00°¢P 





Nan G ee Seaies 


“AW IZW Sc8odil 
a SLO! Vee el ree ere 
YsGdOSde Obs) 
NOCLAGLYLSlG = AQYANS 





OS OP 


SINNDD C421 When 


a2 





oe 6 Oe Sh ee, ae ar 


' 


(eile lasletvalsebiolnsl 


OO°CE 


Ee Whe er aya 


ge eunbt4 


dh 


ee CCl Sie 


WI-WG SEC REL 


LO Ve Ney ees 3) creme co ere ane 


HdgGubSde = Us 
PC ese GU Sis eS SIS raat et 


le Salli bella ee ales) 





SiNNGI GdAaé1 loki 


60 





Pe 


UD “UC 


ge ounbry 


NAW 


US ee 


OOas 0 


Sa 


4 


aS UU" DY 


q34d0Sds Cys 


NOI Nd dls fo 


AD AN 4 





Boas 


| TOWHEN 


a 
i 


Se cls 


6] 





LE aunbt4 


ria 


OD "YG ee OU os Ce UD "PP ae Uee Oe) | OU *0 


; }. 
4 
| iy if 
; 1 | 
of 
i 


Md cb be: [4 CHAN Zia Cie eee 
AoloU INiNaWladdad ——- 04H 


HqedOSds Ubi) 
Eel Tet teats ipsa etal pie 





Bea. PL eh aie Pi Lil Picken! Pel Leticia 


SINNGD O42 7 Tbhish 
62 





OO -To 


oS. 


r 


O 


6 


aN 


r 


r 


ial 


UO °&8 


Vela) | 


OSe ss 


ge aunbl4 


Nal 


00°S8 


J es 


else monet 


NOTINGEYLSTo 


Ob Ge 


@ ee s® 8 +48 


05°08 


» . . + 


“W-HWO PGe ct 
me EL ICT SGI (sll el =) lector clic 
HSGHOSEE Uns! 


mJd4Nn4 


i Pade I ait Bs aia eae 


os 


eee 





SINNDD 04é1 IohiBh 


63 





APPENDIX C - BEAM FOLDING TECHNIQUE 


The theory of Blunck and Westphal [1] assumes that the beam of 
electrons striking the absorber is monoenergetic. The beam of electrons 
produced by the NPS LINAC, or any other linear accelerator, for that 
matter, has a finite energy spread about a maximum or most probable 
energy point. The fact that monoenergetic electrons are not avail- 
able to strike the absorber must be taken into account in computing 
theoretical predictions if a meaningful comparison with experimental 
results is to be made. This has been done here by a technique termed 
“beam folding". 

Beam folding is accomplished by a number of well defined steps. 

Energy distribution curves may be approximated by histograms. These 
histograms consist of a series of "bins" of area W(Q)AQ, where Q is 
the energy loss and W(Q)aQ is the probability of loss between Q and 
Q + AQ. Thus, each bin has an "address", Q, on an energy coordinate 
scale. To accomplish beam folding, the beam distribution must be known. 
This is observed experimentally and approximated in the computer by a 
gaussian curve of appropriate half-width. The energy at which the 
maximum occurs is established by the energy of the beam and the magnitude 
of energy loss incurred by elastic scattering of the beam as it impinges 
on the thin aluminum scattering foil prior to striking the absorber. 
The width AQ is then selected for the predicted distribution. This must 
be a small, but finite number where numerical techniques are to be used. 
This same width is used to break the beam distribution into a histogram. 
This is illustrated in figure 1(a) of this appendix. The reason the 


same width is used for both distributions is to facilitate computer 


programming. 
64 





(a) Beam Distribution 


Counts 





Counts 





Energy 


n Em 
(c) Energy Distribution of E, 


Counts 


EY E, Energy 


(d) Total energy Distribution 





Counts 


Energy 


Figure | 


65 


Each bin of the beam distribution is now treated as a monoenergetic 
beam with energy commensurate with its center location on the «energy 
scale. Each bin also has a definite magnitude, or weight, witl the 
magnitude of the center bin being unity. The formulae of Blunck and 
Westphal is now applied to each of these beams and an absorber distri- 
bution curve results for each, with a maximum amplitude proportional 
to the height of the appropriate beam distribution histogram. Each of 
the absorber distribution curves thus obtained may be thought of as 
being plotted and added to previously determined curves, using the 
energy of the electrons as a correlation point, as depicted in figures 
1 (b), (c), and (d) of this appendix. This is accomplished on the 
computer by adding the contents of each bin of the same address and 
plotting the cumulative total. The bin of address EY is depicted in 
figures 1 (b), (c), and (d). Note that the beam distribution for histo- 
gram E4 is centered over the maximum point of the distribution. This 
is important as a false picture could easily be presented if the histo- 
gram were not symmetrical as the distribution would then appear skewed. 

Since the half-width for a particular absorber increases as the 
thickness (g/cm) increases, the beam folding technique is affected in 
that thick target energy distribution, whose half-width is large, will 
be affected much less by the finite beam distribution than will the 
energy distribution of a thin target whose predicted half-width is not 
more than 2 or 3 times the size of the beam distribution half-width. 
Thus, in selecting bin width (aQ), it is desirable to make the width 
proportional to the target thickness so as to apply the beam folding 
technique in greater detail to thin as opposed to thick targets. The 


computer program (Appendix D) applies the target thickness in direct 
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proportion to establish a bin width for subsequent beam folding 
calculations. Thus AQ is much smaller for the thinner absorbers 


compared to thick ones. 
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APPENDIX D. Computer Program 
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TSTPTRUTIIN OF BFAM ENSOGY 
PASE DON THe EXPERIMENTALLY 
ee Dee ee ey eae) Eee (RE AM 
GEAR ee tn eS Ss 


FUNC 


a> 
D 
on 


n 
a WN 
ZZII=ZVN=E NA vA NS «KK ODNAN TM NA BMNIAN™ 
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Thame 
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‘ereGe na” AS A "BEAMN, THE ENLEORING PORTION OF 
feaeuieS AN ENERGY MISTRIBUTTION FOR EACH 
FAMS" AND STTVRES THEY IN "RINS" et aa am 
fee MAT TVS AME UN SES eRe ere aN COUNTS IN 
WOT, 9 TS THE FNE®GY LOSS OVER WHICH THE R=w 
IS 


INTEGRATED ANY) NX US eM AY WATCH 
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TEGP ATION, et SGN 
. 


TS CAPRIFN OUT IN SUBROUT 
x= ae 


CEO ion seh? SEACH SUICCESSTVE 
INF 


eMC) 


Cooma? 171771 35.)* *5 50=NX 
re( a, L*. ET IG) may nos 
N=9 
GOTOLONS5 
Bees [F(QsLT.%)1GNTNIN0G 
FIT=F (1) 
Migea o 20ODC}~O 
VALIM=.75%9 
PN NOG37(05,%L Se oe legless 2) 
Sammie es (XL, VALIM, ECT. 78.51.0057, Ay?) 
beac s7 (VAL IM,O »FCT, ZC,E1.0.T. As Z} 
WQ=7A4+7B47C 
WOT(LEN=-T) =WOT (EL tN=-1 D4 WOFNE (IL) 
XX(L4IN-1T =F (LI-0 
; v= Et N=] 
C THE NISTRIBRUTION FOR THE HISTOGRAM CORRESPONDING TN THE 
c SeNTER OF THE BEAM DISTPIBUTION IS REMAVED AND STIREN IT 
. WUF ANI XUF AS THE UNFOLNEDN COUNTS AND ENERGY 
: RESPECTIVELY, 
eleewN Gs TCENTIGOTOI9NO5 
WHF(N) =WO 
XUFUNI=FI-9 
K=N 
GOTO 1°95 
fO06 CONTINIEF 
. Oe sig y=) ,7 
Seer NATE tOW NUMBERED *"RINS"” AS THEY ARE ONLY PARTLY FIL= 
a LED AND ARE NAT REPRESFNTATIVE JF THE CUMULATIVE 
e DISTRIBUTION, 
WUFC(C T) =WHF (8) 
1717 WOT(T)=WOT CR) 
DY 6 T=1,4 
6 TI(TI=TSA(T) 
GATO 779 
; 177 00778 Ne=l,K 
c pe TS 3 THE SHINENLNFEN THEORY CURVE. SeeetNG EXX( 2 ) SOUAL 
cc ZEPM ALLOWS PLOTTING OF THIS THENRY CURVE WITHOIIT ANY 
EXPERIMENT AL DAT A, 
WOT(N) =WUECN) 
XK (NJ=XUFON) 
KSTG=1 
J=K 
778 CONTINUE 
EXX( 1) =9 
NI 5 T=1,4 
Ser stl y=1 ol} 
779 2 ae 
‘ DONICOT M=3, 
C METERMINE THE MAXEMUM POINT NE THE OTSTRIRUTIGN, WOTMX, 
c AND THUS SSTARLISH THE YOAST PROBABLE FNERGY LOSS. QP, 
§ 
se WOT (N=1)5 AND, WOTON-1),G6EF,.WOT(CN=-2)) 
1G9TNI1ON 
GATN ae 
1NCB FECWOTINALILELT,WOTMXIGOTA LOC? 


Brana 1 
OP=FI=xxX(N~?t) 
mee? SCONT INUE 


C 
C NORMALIZE ALL POINTS OF OTSTRIRBUTTON AND STORE ITN HW, 
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N9e1 noc Me]. J 
HW(MJ=HWOT OM) /SWOTMYX 
INCSO CONTINUE 
jay—) 
HWI=0 
C DFTFERMINE HALFWIDTH OF THE THEARETICAL NISTRIRBUTIAN, HWA, 
No TATA ova tlew 
TPC HWOM41) .GTe,5,AND,HW(M).LIT. 5) HWITH=XX(M) 40,56- 
THW(MY) ENELX/(CHWOMET)-HWOM) } 
ITF CMSFQ, Jo ANDZHWI, 5959 IGNTI1IIS 
Gone i271 
11° WRITE(46,?9) 
HWA=1 
Gone Ten : 
20 FORMAT(/,2X,*ND HALEWIDTH OBTAINED ?Y) 
1721 TE CHWOM).G7,55,4ND,HWOM41)5LT, 5 SPHWARXY(M) CHW (MY) 
Y-, S)V*NEL XS (HWOM) -HW( M41) )-HWT 
1910 CONT INE 
: Bee ISH CUTDEF STGNAL FOR PLOTTING SUBROUTINE, 
; RE J+) y= 
: ESTABLISH CORRFCT GRAPH TITLES FOR ABSNPRBFER USED, 
[F(7> Greate GMa: o 
DN (Rea 
3 Poti )er att) 
G@T0101 1 
Poe Te( 266! 12951)G9TN50 
DN 1MI=1,4 
memeee 1 j= TOC] ) 
GONTOIOLI 
BQ 1F(72.67T.503,;1)GAT 64 
Hea 8 1=3, 4 
aw (1 }=TS (7) 
Game ol? 
Cote (7. G6 eho )GOTNS 2 
NO 9 T=1,4 
9 T2(T)=T9 (T) 
G@rQ1o0)1 
Aon 7 15144 
: Peeuet tl jy=atT? ty) 
e SEND DATA TON PLOATTING SUBRNUTINE FOR PROCESSING, 


merry CALL GRAPH (XX, 
TCORR ,CTSNAR, OP) 


WOT EXXeEXVe TS, T 32, TA oT sHWASZ FI a Ty 


HALF ee. 


C 
C PRINT RESULTS OF COMPUTATIONS IN TARULA® FORM, 
c 

WRITE(6,722) OP,HWA,7T,7,F I 

TF(KSIG.=0Q,°) GATA TIT 

TUE ORMAT (7 2104.0 ) 
CPR Me ext OPS § F754, * 
1*®* TCGM/S9O-C4) 2 §,F65,%,° 75 
296 O68 25/3 

STOP 

END 
C 
C 
C 
c 
¢ 
C 
C 
C 
C 
s 
C 
c 
C 
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io GIG GC? CIGIte NOOO eo 


OAIMAAOASY 


Ge 


SAL ElINCTION WEON SCX ¢51ete79A,9° ) 
eases SUNG TON CALCULATES THAT PART OF THE FENFPGY 
Maison TON Te AN ELECTRON DUF TT LASS AF ENERGY BY 
IONIZATION WHILE PASSING THRRUGH AN ABSNQRER, 
Re eee ie SAS 7) 


Meee THE AVERAGE TONTZATTOR PXTENTIALS PER FLECTRON FPR 
er) FOR FACH SHELL. BEGINNING WITH FHF K SHELL, 


oreo HE NUMBER OF INNI7ZATION POTENTIALS USED FOP LEAD 


etrtS CONTEXT, L=LFAD, G=GANOLINJUM, S=TIN, C=CNOPFR, 
AND A=ALUMINU", | 


REAL*R [L06)/5 9880090 15 914709, ,0CR9NN,, INF494, 
ee OO 754, OFC IO? / 
RFAL*& Ney ee egy se 77g 1 as 4 
Bees SG oy pnt 4 5 99795 44 10 Bh, 5900 24,,0909247 
Bee 8 N60 9 7 2 Bo gp ERS 47 55 
REAL*B TS(E)/,9292,,90444,900774,099162, 900001 / 
emer NN St DPS 25 98s pl 8a ela g 
Were A 1C( 7975 00898,,0°0069, ,.090%68/ 
REAL *S NC(3)/27,,8.,195/ 
ponies a2 175 O96, ,OCORA/ 
RFAL*8 NA(2)/2,,8, / 
NPOTL=6 
NPATG=5 
NPOTS=5 
NPOTC=3 
NPT A=2 
Soo SUMMATION OF ITONIZATION POTENTIALS FOP THE 
EXPERIMENTAL ATOM, 
SB=9 


XMASE IS THE REST MASS OF AN ELECTRON IN MEV, 

XMASE=, 511906 
RETA 1S THE NORMAL V/C USFD IN RELATIVISTIC CALCULATIONS, 
GAMMA IS THE NORMAL TERM USED IN RELATIVISTIC CALCULATIONS 


BETA=DSORT(CET*EI+ 25% 19 XMASFI/ (ET EL +2, STS XMASE 
1+XMASF XX4ASE) ) 
GAMMA=1,/(15-BETA#*2)**(1,/25) 


AT THIS PNINT, THE PRAGRAM MUST COMPUTE THE DRAPER CONE] 
STANTS FOR THE FXPERTYENTAL 7, A SEQTFS NE LOGIC STATF- 
MENTS SELECT THE CORRECT ENARMULAE, 


13 IF(7.GT.13,1)IGOTN?9 


Seve ote e SEXPERIMENTAL Z;5 SS eer AGT ON POTe N= 
Tae ONES SHEL: 6 


TA 
V=ALAG (73 “E1/ CIRC (14-8 e rae?) ) 


Pee een ee iat )F= 


2B? AND TALS THES 
TE THE AVEPAGS FNERGY 
rt 


EN H 
TS UST) TO COMPUTE 
ION AND ARE SIHPPLTE 
PE RB AND Ay RESPECT 


BV THE STEQNHEI[ MES 
VELY, F2NM HIS PAPES, 
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NMOAOAOOADVEY 


MOANA AI) OO) 


OO 


AVEPAGE FNERGY LASS DUE TO IANT7Z ATION, NUMBERS HEOEFTA 
REPRESENT OTHER STEPNHEYMER CONS TANTS LS Se Oe WSs 
Seaek A «M996 
SMALL ™M ae 
Xx] Aa 
-C oT 
DEL T =45 AN 6*ALAGIOC( RE TAX GAMMA )-4,214,N9I06K* 
103, -DLOGIA(CRBRETA*GAMMA) )«*2,5]1 
Gen). 2 11 
Poteet lel ls 279R)YGNTNS49 
Nn 12 l=]. Nea C 
RB=IC (LI*NC(T)#0LOG(2, MEI/CICCIY*(1.-RE TAs ©?) ) } 
72 SR=SR48 
B2=15,9 
ae) 
Ost —-4, 696 OU NGIO( SETA -=GAMMA)-4,744+,11908 
10346 -NDLAGITO(RETAXGAMMA ) )* «3,28 
Comm 17 
S0 ier a>) os O05) 1GOTOGRG 
DN 732 Y=l,NPOTS 
REIS(CTITENSCTI*YOLOG( 2,7 FI/CISCTY*K(1,-BETA#*7))) 
73 SR=SKR+BR 
B=] 3, 33 
A=, 064 7 
NELT =4,69N5*“NLOQGLO(CRBETA*GAYMA )-6, 28 4,494%(2,- 
LOLOGIO(BETA* GAMMA) ) 
1**2, 52 
Ga 1 311 
64 1TF(75675,64,1)GNTNAR? 
ON 74 eS 
B=IG(T)Y*NGCLTI*OLOG(7,* FT/SCIG(I)#(1,-BETA**2)) ) 
7T4 SR=SPHER 
R2=13.23 
Al=_,96274 
NELT=4., 4694 *NLOGIAC(BETA*GAMMA )-5, 94,41 8"(2,- 
LDLOGIO(CBETA*GAMMA) ) 
1*#*?7?,19 
Caimi 1) 1 
Bre ele Ol wt 5 825) )GOTOR ON 
NO 75 T=1,NPOTL 
B=H=ILCTY¥NELCTP*OLNG(2,* EP/CILCT*01,-BETA**2))) 
75 SB=SBEB 
Boose! 
Al=,06%8 
DELT=4, 6NG*ENLOGIAN( RETA* GAMMA )-6,934,95527% 
1(4.-DLOGIOCRE TA? GAMMA) )¥*2,41 
GO TO 1l1t 
AR TS THE OUANTITY SMALL A, FROM BEUNCK ANN WESTPHAL, 
MULTIPLIED BY R, THE TARGFT THICKNESS IN CENTIMETERS, 


1TJ]1 AR=9,154*7*T/A/BETA/BETA 


PeMiPUTe BS 2, THE SMALL PB SQUARED FROM THE SLUNCK AND 
belS =GANG PAPER. 


BS 2=36 “S57 (/7*A"®) 


VAR1T, VAR?, AND VAR ARE CANVENTENCE STORAGE LOCATIANS 
WSED TH STORE SEM OSCR ives aes ert Meet SAV Ce AGG 


ENERGY LOSS, 
VARTSAL*T/ BETA 2 
VARP=BR744,4 247, *NLNAG (RE TA*GAMMA ) 
VAR3=NLIG(ET)—-RETA*E? 
Save S THeeeI NAL VALUE FMR THE AVES AGE ENFP GY LOSS. 
QAVE=VAR1L* (VARZ4EVARI-DFEL® ) 
BL=DSQRT(BS?) 
Meeseete PORTION MF SAME TOTAL ENERGY LOSS QO NUE T9 
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Crt CC ae 


AAAAD 


REAL FUNCTION WRAN7A(XK,F Ty T7540) 
ements UNG Tet CALCULATES THAT PART Af THE EMERGY 
ee aN iaeN ELECTRON DUE TS! LSS MF ENERGY RY 
RANTATION WHILF PASSTNG TH8TUSH AN ABSORBER, 
IMPLICIT REAL*R(A-HyN~7)} 


U GA Peer STORAGE LOCATION USED IN COMPUTING 
ALPHA, 


=e 5 O/7*~ (16/450) 
ear VALUE OF ALPHATS FOR TARGET THICKNESS, ATOMIC 
NUMBER AND ATOMIC WEIGHT OF EXPERIMENTAL ARSMPRER, PB IS 
ince Ge | THICKNESS FXPRESSED IN CENTIMETERS, 
ie te OO) ee 7 PAY ( 4. OO OR OLCG(UT 4198/9, 7) 
SOMPUTE VALUF FOR BR, A NOPMALIZING CONSTANT, 
B=1. /DGAMMA(ALPHP #1, ) 
meee SWRA, TRE MNUMRER OF COUNTS =XPECTEN AS A PESIIL® OF 
Seem een Se GY ETSSES, X BEPRESENTS FHAT PARTICN OF 
eer ENERGY (SS 0 WHICH [1S NOT LOST By TONIZATTION 


WRAD=(B*ALPHRE (X**ALPH®)/(ET** ALPHRII/X 


RETURN 

END 

REAL FUNCTION BEAM*R(FI,HW,X) 
IMPLICIT REAL*P(A-H,9-Z2) 

a 

ALFEA=4,*DLOGCY )/(HWHW) 
BEAM=DN=XP(-(CX-FI)*%47) “ALEA) 
RETURN 

END 


ROAL FIINCTICN FOT*2(X% 7 sANg TVET OO) 


PURPNSEs THIS FUNCTION MULTIPLICATTVFIY JOINS TOGETHER THE 
PRENICTION NF ENERGY LOSS NYE T2 TONTZATION ANN THE L9SS 
NUE TO RADIATION, 


PLICTT REAL*¥S(A-H,9-7) 
)* 


T=WRAD(X FI i are 


] 

Ee WION(X SFT + Te+725A09 ) 
RFTUPN 

EN 


N 
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